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Abstract 

We study the Generalized Brans-Dicke cosmology in the presence of matter and dark energy. 
Of particular interest for a constant Brans-Dicke parameter, the de Sitter space has also been 
investigated. 
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1. INTRODUCTION 



The Brans-Dicke (BD) theory is defined by a scalar field and a constant coupling 
function u as perhaps the most natural extension of general theory of relativity that is 
obtained in the limit of a; — oo and ip = constant. The theory appears naturally in 
supergravity theory, Kaluza-Klein theories and in all the known effective string actions. In 
FRW cosmology, the theory gives simple expanding solutions for scale factor a{t) and scalar 
field (p{t) which are compatible with solar system experiments |l| 

The generalized BD theory, sometimes referred to as graviton-dilaton or scalar -tensor 
theory, is instead, defined by u which is implicitly a function of time u){t). Naturally, a 
few attempts have been taken to study the dynamics of the universe using this formalism 
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e during the early 
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In this case 



The belief that modified gravity theories may have played a crucial ro 
universe has recently been renewed by extended infiation (for example see 
a scalar tensor gravity theory allows the first order phase transition of the old infiationary 
model to complete. This arises because the scalar field that is essentially the inverse 
of the Newtons gravitational constant, damps the rate of expansion and, in the original 
extended inflationary model based on the BD theory, turns the exponential expansion found 
in general relativity into power law inflation [3]. However, BD theory is unable to meet the 
simultaneous and distinct requirements placed by the postNewtonian solar system tests and 
by the need to keep the sizes of the bubbles nucleated during inflation within the limits 

n n 

permitted by the anisotropics of the microwave background [8| [9j. 

In order to carry out a detailed study of the dynamics of the cosmic evolution in this 
formalism, knowledge about exact time-dependence of a(t), ip{t) and uj{t) and energy mo- 
mentum distribution in spacetime is essential. In part of this paper, similar to lOj, we have 
obtained red shift-dependence of uj{z) with the power of red shift z determined in terms of 
the exponent of scale factor a(t) which is taken to vary as a(t) oc and equation of state 
parameter for the matter and dark energy contribution. With the help of observational 
evidence we obtain certain information about the parameters describing the cosmological 
model in particular regarding the early and late time behavior of the universe. We also in- 
vestigate both empty and flUed de Sitter space case with constant BD parameter and shows 
that the result is consistent with recent measurements. 
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2. THE MODEL 



We consider a flat Universe fllled with pressureless matter and dark energy both described 
by perfect fluid. The fleld equations in generahzed BD theory with time dependent u, are 

17-2 2 I IT- • ^ -2 + Pm 

H If + Hifiif -Qf= ^ (1) 

2Hip'^ + SH'^ip^ + + "iHipip + i:pip = -Px<^, (2) 

where px = axPx, Pm = ctmPm are the equations of state for dark energy and matter and 
the scale factor and scalar fleld are respectively a{t) and if{t). In addition, the equation of 
motion for BD scalar fleld is given by 

.. , ^rj. Px + pm- ^Px UJ^ 

From equations ([1]), (j2]) and ([3]), the energy conservation equation can be obtained as 

{Px + Pm) + 3H{px + Pm + Px) = 0- (4) 

Note that the wave equation for the BD scalar fleld, ([3]), is not an independent expression 
as it follows from the Bianchi identities alongside equations ([T]) and ([2]). In addition, the 
dynamics of the scale factor is governed not only by the matter and dark energy, but also 
by the BD scalar fleld, (p{t). 

One may assume that the matter and dark energy interact with each other, thus the 
growth of one is at the expense of the other. Then the conservation equations for them are 

Pm + ^Hpm = Q, (5) 
Px + 3H{1 + ax)px = -Q, (6) 

where Q > stands for the interaction term. Alternately, one could construct the equivalent 
uncoupled model described by: 

Pra + 3H{1 - am,eff)Pm = 0, (7) 

Px + 3H{1 + ax,eff)Px = 0, (8) 

where 

orlpm 

Q / ^ 

ax,eff = ax + ^7T— • (10) 
6Hpx 
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The wave equation ([3]) is not altered by the interaction equations ([5]) and ([6]), since although 
the matter and dark energy components do not conserve separately the overall fluid -matter 
plus dark energy- does. Thus, one may introduce the total energy density ptot = Pm + Px, 
and from equation ([7]) and ([H]), obtains 



Ptot + 3if(l + atot)ptot = 0, (11) 



with the solution 



where 



Ptot cx a-3(^+'^-'), (12) 



<ytot = — — — = c^x^x, (13) 

Pm + Px 



and Qt = —. One can also find the rate of Qr as 



^ _ -^H{a^^eff + am,eff)PxPm _ SHatotPm + Q ^^^^ 
Ptot Ptot 

or in terms of red shift z, 

^, Hax,eff + am,eff)Pxpm _s 

= AJ^) • '''' 

where ^ = 1 + z, and a = 1 is the present value of the scale factor and ""' means derivative 
with respect to z. 

Also in terms of the red shift z, equation ([8]) can be rewritten as 

, ^ 3(1 + ax,eff)px , . 

Px , 

where the equation based on the sign of 1 + dx^eff shows whether the density of dark energy 
will increase or not as the red shift becomes low. For positive sign, the density decreases 
like the quintessence, for negative sign, it increases like the phantom, and when it is zero 
the density is invariant like the cosmo logical constant. 

3. THE GENERAL u 



Following paper |2| for a oc t^, if oct^ and time dependent u, equation ([T]) gives, 

P 
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One can rewrite equation (fT7|) in terms of red shift z, 



where its derivative with respect to z is given by 



After some calculations, we can also rewrite equations ([I]) and ([3]) as 

,a ^ 2 _ (2w_+3)vj^ 



ilT + ~) I'^f = Ptot, (20j 



/3[^ii^^/3 + ^^i^^] = 0. (21) 

From equation (I2T!) one finds that for 7^ 1, /3 restricted to be or —25. In case of 
atot = 1, there is no constraint on /3. For /3 = 0, from equation f|T8|) one finds that a; — ?■ 00 
and from equation (!20|) we obtain ip = ipo = constant and a (x t^ . So for /3 = 0, Brans- 



Dicke model goes over to General Relativity 



11| and to obtain 5, one has to solve equations 



of General Relativity [12|. In case of /3 = —25, equations (fTSll and ( fT9l) reduce to. 



a;(^) = _-I^(l + z)a+3"-)-2/^, (22) 
2o^ 

It is clear from equations fl22l) that, the parameter atoj which takes different values in different 
era, controls the z dependence of u in different era. For today value of 2; = 0, we have 

-0 = ~. (24) 



and 



, -^(l + 3a,o,) + 2 

^0 = 1^, , (25) 



where for the present acceleration of the universe that 5 needs to be greater than one, ujq 
given by f l24|) has the minimum value of 1/2 in agreement with the observation [isl. 

Further, in the last scattering surface, during the galaxy formation era (1 < ^ < 3) 
where dark energy density must be sub-dominant to matter density {atot > —0.5), we have 
-1/12 <cu < -1/14. 

In the Big Bang Nucleosynthesis (BBN) era where the presence of dark energy should 
not disturb the observed Helium abundance in the universe {{atot)BBN > —0.21 at 2; = 10^°) 



|l4l |. we have —9/32 < u < —9/128. This also shows that at sometimes in the future, 
z = —1, we have a big rip and u — > — oo. 
One also finds from 022]) and that 

- = {{l + 3atot)-2/6){l + z)-\ (26) 
oo 

where the ration for today is negative, in the distance future for z = —1 it goes to minus 
infinity and in the distance past where z — )■ oo, it approaches zero. 



4. DE SITTER SPACE TIME WITH CONSTANT u 

We now assume an empty de Sitter spacetime with the solution H = Hq. Then, for 
V? oc and constant BD parameter u, equations ([1]) and ([2]) can be solved: 

H', + Hof3 - = 0, (27) 

3H^o + 2W + (| + = 0, (28) 

to give 



(-2 ± v/-8 - 6uj)Ho 

^ ' ^^^^ 



P ^ (3±v^)^o. (30) 

CO 

For these two solutions to be consistent implies that oo = —4/3 or u; = —3/2. For equations 
([1]), (|2]) and ([3]) to be simultaneously satisfied only oo = —4/3 and so (3 = —3Hq is acceptable. 
For large value of Hq, during infiation, while the universe expands exponentially, the BD 
scalar field drops exponentially. 

In the presence of matter and dark energy we may also have a de Sitter solution H = Hq. 
Then the equations ([1]) and ([2]) become 

+ H^f3 _ ^f3^ = ig[-3Ho(i+mo*)-/3]i^ (31) 
6 3 

3H^ + 2HoP + (^ + l)/32 = -a,„,e[-3^°(i+"-')-'^l*- (32) 

These equations are satisfied when 

H',+Hof3-^f3' = ^, (33) 
+ 2HoP + (| + l)/3' = -atot, (34) 
-3Ho{l + atot)= (35) 



Using equation (l35l) in (133!) one gets 



Similarly, equation ( 135|) in conjunction with equation (!34|) gives 

_ 2(3gg(l + 2atot) - mlil + fttot)^ - atot) 

For equations ( 133|) . and (!35|) to be simultaneously satisfied, above two values of u should 
be equal. Imposition of this condition leads to 



iSHi + 1) ± ^9H^ + A2HI + 1 

From the above solution we find that for negative sign and Hq > 10, atot = —0.33 or for 
Hq < 0.15, atot > —0.33. For positive sign, for Hq ^ 1 or Hq 0, we have atot — 0. In 
case of atot > —0.33 or from equation f fT3l) . < 0.33, this is consistent with last Scattering 
Surface, during the galaxy formation era (1 < 2; < 3) where dark energy density must be 
sub-dominant to matter density and accordingly < 0.5 . Then one gets u = —1.5 and 
13 = -142.7. 

From the above argument as to the Brans-Dicke parameters u are concerned, it is negative 
and of the order of unity. This could be considered as an unsatisfactory result, in view of 
the high lower limits imposed to u by astronomical tests in the Solar System. 

A possible solution of this contradiction as discussed is in considering a non-constant 
coupling function u{t) in Generalised Brans-Dicke theory. Thus, the value of such a function 
can change with the cosmic time and, in the limit t — )■ 00, it could agree with local measured 



values 



15| . This argument is based on the scalar-tensor theories in which u depends on the 



scale, being very high in the weak field approximation of Solar System that probe only a 
limited range of space and time. To effectively constrain more general scalar-tensor theories, 
one would also like to have strong-field experiments, such as that provided by the binary 



pulsar 16|. It was also pointed out that in cosmological models based on more general 



scalar-tensor theories in which u can vary, there is generally an attractor mechanism that 
drives w to 00 at late times |17 |. 
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5. CONCLUSION 



In this work,we have derived the exphcit time and red shift dependence of the Brans- 
Dicke parameter uj by solving gravitational field and wave equations of generalized BD theory 
consistently in the present of matter and dark energy, assuming power law behavior for the 
scale factor a{t) and scalar field ^{t). Similar to the work done in j^, we find two consistent 
solutions of the field and wave equations. One solution leads to General Relativity and the 
other one leads to a ^-dependent uj{z) whose red shift dependence is governed by the equation 
of state parameter atot- Consequently, oj{z) exhibits different behavior in different epochs of 
the evolving Universe characterized by its dominant matter/dark energy components. We 
also find that the ratio ^ is a negative monotonically increasing function of z. In particular, 
for an expanding universe, we have studied the empty de Sitter space with constant u and 
find that in the era of inflation that Hq is high, the scalar field drops exponentially and 
u = —4/3. Moreover, in the presence of matter and dark energy, we also find that u is 
negative and we are able to explain the last scattering surface constraint 
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